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We theoretically investigate the influence of designed pulse sequences in restoring quantum co- 
herence lost due to background noise in superconducting qubits. We consider both 1/f noise and 
Random Telegraph Noise, and show that the qubit coherence time can be substantially enhanced 
by carefully engineered pulse sequences. Conversely, the time dependence of qubit coherence under 
external pulse sequences could be used as a spectroscopic tool for extracting the noise mechanisms 
in superconducting qubits, i.e. by using Uhrig's pulse sequence one can obtain information about 
moments of the spectral density of noise. We also study the effect of pulse sequences on the evolution 
of the qubit affected by a strongly coupled fiuctuator, and show that the non-Gaussian features in 
decoherence are suppressed by the application of pulses. 



I. INTRODUCTION. 

Quantum decoherence, the continuous decay ("the loss 
of quantum memory") of a quantum state due to its in- 
teraction with the environment provides the conceptual 
connection between the microscopic quantum and macro- 
scopic classical worlds. 1 ' Understanding and preventing 
decoherence is also central to the fledging field of quan- 
tum computation, as the loss of quantum coherence leads 
to errors in the processing of quantum information. In 
fact, quantum error correction protocols, necessary for 
quantum computation, require the decoherence to be be- 
low a certain threshold. 2 It is therefore of paramount 
importance that the decoherence of qubits, two-level sys- 
tems used to store and process quantum information, 
is suppressed as much as possible. In this article, we 
develop realistic strategies, based on the application of 
designed external pulse sequences, that suppress an im- 
portant source of decoherence in an imp orta nt class of 
qubits, solid state superconducting qubitsP^ The deco- 
herence mechanism considered in this work is that due to 
classical noise, i.e. a situation where the qubit couples 
to a random classical temporally fluctuating field. Such 
noise is, in fact, the major s ource of qua ntum dephas- 
ing in superconducting qubitsP 6 * 7 * 8 * 9 * 10 ' 1 -^ and therefore 
the pulse sequences proposed in this work should be use- 
ful in restoring coherence in solid state superconducting 
quantum computer architectures. 

For many decades, in the field of magnetic resonance, 
pulse sequence techniques have been s tudie d as a method 
of reducing spin ensemble dephasing ! 12 * 13 ' The most fa- 
mous sequences are Hahn's spin echo, Carr-Purcell- 
Meiboom-Gill (CPMG) sequence,^ and periodic dynam- 
ical decoupling (PDD). Spin echo (SE) is the simplest, 
consisting of just one it pulse, whereas CPMG is its multi- 
pulse generalization. The PDD sequence, introduced in 
the context of quantum computation, was designed to 
average out the influ ence of the environment, effectively 
decoupling the qubit IHUSUSi Dynamical decoupling was 
developed further by introducing an idea of recursively 



defined sequences,^ termed concatenated dynamical de- 
coupling (CDD). Recently, using a spin-boson model 
Uhrig obtained a new sequence-^ (termed here UDD), 
which nearly completely suppresses short time decoher- 
ence under certain conditions. UDD was later shown to 
be universal,^ in the sense that for for any pure dephas- 
ing Hamiltonian the n-pulse UDD sequence leads to the 
cancellation of n orders of the time expansion of the off- 
diagonal element of the qubit density matrix p + _ (t). 

In this paper, we study the effect of pulse sequences 
on the decoherence in superconducting (SC) qubits, sub- 
ject to classical 1/f and Random Telegraph Noise (RTN). 
Experimental studie d 5 * 20 ' 21 * 22 ' have shown that SC qubits 
suffer decoherence from 1/f noise, which is associated 
with fluctuations of electric or magnetic dipoles in the 
insulating materials. In charge qubits, where the area 
of the tunnel junctions is small, it has been established 
that the qubit is often coupled to a few two-level fluctua- 
tors (TLFs) , which can be treated as classical sources of 
RTN.' 5 * 6 ' 9 ' 11 ' Here we focus on the case of charge qubits, 
in which the charge noise is dominant, but our results, 
with minor modifications, are applicable to phase and 
flux qubits. 

It has been experimentally shown that the coherence 
of the Cooper-pair box charge qubit is significantly pro- 
longed by the application of the SE sequence.' 5 ' 11 * 23 ' Char- 
acteristic plateaus seen in the echo signal have been ex- 
plained theoretically as arising when the noise is domi- 
nated by a single clas sical TLF coup led to the qubit P 
Beyond SE, only PDD 1 24 ! 25 ^ ! 27 ™ and CPMCpH se - 
quences have received theoretical attention in the context 
of SC qubits. In this work, we suggest the use of more 
sophisticated pulse sequences, such as CDD and UDD, to 
suppress noise-induced decoherence in SC qubits, finding 
that, depending on the details of the noise, CPMG or 
UDD is optimally effective in reducing decoherence in su- 
perconducting circuits. We emphasize that earlier work 
in the literature on CPMG, 1 ^ CDD' 31 ' 32 ' 33 ] and most re- 
cently UDD 19 pulse sequences was carried out entirely in 
the context of electron spin decoherence in a nuclear spin 
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bath, in which the quantum correlations within the bath 
are important. On the other hand, the bath fluctuations 
due to charge noise in superconducting qubits can often 
be treated classically as established in Ref. |5] 

Apart from prolonging the coherence time, pulse se- 
quences could be used to gain valuable information about 
the environmental noise, as the time dependence of de- 
coherence is different for various sequences. The micro- 
scopic origin of the noise affecting the coherence and en- 
ergy relaxat ion in the SC q ubits is still a subject of ongo- 
ing research j 34 l 35 l 36 l 37 l 38 l 39 l anc | we discuss here how pulse 
sequences can be used to learn more about the noise spec- 
trum at low frequencies. This approach is complemen- 
tary to using the measurement of energy relaxation time 
of the qubit for noise spectroscopy at high er frequencies, 
of the order of qubit energy splitting? 40 * 41 ! 

We consider here the experimentally relevant ^ * 20 * 21 ! sit- 
uation in which decoherence is dominated by pure de- 
phasing (i.e. "T2") processes, and not by energy relax- 
ation (i.e. "Ti") processes, Ti^>T 2 . In the current ex- 
periments pure dephasing is dominant mechanism of de- 
coherence away from the so-called optimal bias point fSBU 
The decoherence at the optimal point in present charge 
and flux qubit designs is limited by T\ processes (i.e. 
T2 ~ 2Ti ) , since the effect of noise is then suppressed to 
the first orderP However, the requirement of keeping the 
qubit at the optimal points at all times might be overly 
confining for a system of multiple interacting qubits. 43 
Thus, the ability to prolong the coherence of the qubit 
in the pure dephasi ng regim e is still desirable. Further- 
more, in phase qubit d 24 ' 44 ' 45 ' there is no optimal point and 
the effect of pulse sequences should lead to a substantial 
increase of T2 in the c ase wh en qubit coherence is T2 lim- 
ited. Another strategy^ 46 * 47 ! for suppressing the influence 
of the noise on the qubit was implemented recently.^ In 
such a "transmon" qubit the coupling to the charge noise 
is exponentially suppressed. However, this qubit is still 
sensitive to the flux noise, the relevance of which becomes 
prominent away from the optimal flux bias point. 

The article is organized in the following way. In Sec. [H] 
we introduce the pure dephasing Hamiltonian and de- 
scribe the types of noise which we shall consider. Sec. |III| 
contains the overview of various pulse sequences appli- 
cable to the pure dephasing case. In Sec.[lV]we present 
the analytical solution for decoherence under pulses for 
the case of Gaussian noise, and we discuss how the pulse 
sequences act as filters suppressing the influence of low- 
frequency noise on the qubit dynamics. The calculations 



for Gaussian l/f a noise are presented in Section IV B 



where the role of the ultra-violet cutoff in the noise spec- 
trum is highlighted. We also introduce the idea of using 
the UDD sequence to obtain the quantitative information 
about the low-frequency noise spectral density. Finally, 
in Sec. [V] we present the results for decoherence due to 
the RTN. We identify the regime in which the application 
of even a few pulses leads to the increase in the coherence 
time, and we find that the analytical Gaussian approx- 
imation to calculation of decoherence is asymptotically 



exact for large number of applied pulses. 



II. THE HAMILTONIAN AND THE MODEL OF 
THE NOISE 

The limit of energy relaxation time being much longer 
than the dephasing time, T\ 3> T2 , corresponds to using 
the pure dephasing Hamiltonian to describe the qubit- 
environmcnt interaction: 



H=i[n + /3(t)]*, 



(1) 



where f2 and j3(t) are, respectively, the qubit energy split- 
ting and a classical random variable representing fluc- 
tuation of the energy splitting due to coupling to one 
or many TLFs. The function j3(t) represents a classical 
stochastic process, given by 



(2) 



where £j(t)=±l/2 corresponds to the RTN signaP^from 
the i-th TLF, with m being the corresponding coupling 
strength. 

The stochastic processes are defined by their corre- 
lation functions. The two-point correlation function is 
given by 



S(t 1 -t 2 ) = (f3(t 1 )0(t 2 )) 



(3) 



where (...) is the average with respect to the noise re- 
alizations, and we have assumed here (/3(t))=0. The 
Fourier transform of the two-point correlation function 
is the spectral density of noise (more generally referred 
to as the first spectral density, see Ref. I5U]| : 



S(u) = 



(4) 



When the statistics of fluctuations are Gaussian, the 
noise is completely defined by S(t) , and the average over 
the noise realizations can be written as a Gaussian func- 
tional integral over all possible realizations of j3(t) : 



= V(3e 



-i/dt 1 /dt 2 /3(t 1 )S- 1 (ti-t 2 )/3(t 2 ) 



(5) 

with S* -1 is defined by 

dif'S' 1 ^ - t")S(t" - t') = S(t - t') . (6) 



On the other hand, when the noise is non-Gaussian, one 
has to consider also higher order correlation functions. 
For the relevant here case of the RTN, the two-point cor- 
relation function and its Fourier transform are given by 



S(h - 1 2 ) 



^1)^2)) = 

2 

7^ 



4 7 2 
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(8) 
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where 7 is the rate of switching between the two values 
of £ = ±1/2. We have used here the high temperature 
limit {k^T ^uS) for the spectral function, i.e. the noise 
is symmetric] 9 ! 49 ! with both rates of transitions between 
the two states of the TLF being equal to 7. If the noise 
had Gaussian statistics, one would have been able to ex- 
press the higher order correlators through the two-point 
correlation function S(ti — £2)- This is not true for the 
RTN, and we refer interested reader to Appendix [A] for 
more details. 

When many TLFs with a log-uniform distribution of 7 
(i.e. with probability of finding a fluctuator with a given 
7 being P(j)oc 1/7) contribute to (3{t), the spectral den- 
sity is S(u))<x 1/lo. This is the well-known 1// noiseP^It 
extends to an infrared cutoff frequency Wj r below which 
the spectrum flattens out, with values of w; r /27r < 1 Hz 
in SC qubitsP^ The log-normal distribution of switch- 
ing rates arises when 7 depends exponentially on another 
quantity having a uniform distribution. For example, in a 
model of localized TLFs one obtains a log-uniform distri- 
bution of tunnel splittings which depend exponentially on 
the tunnel barrier height,^ and in th e rece ntly proposed 
model of the Andreev fluctuator batrPSEU th e switching 
rate of the effective TLF depends exponentially on the 
distance between the pair of impurity sites participating 
in the Andreev tunneling process. 



III. DECOHERENCE UNDER PULSE 
SEQUENCES. 

We consider the decoherence of the qubit prepared 
initially in the coherent superposition of its "up" and 
"down" states. Specifically, at time if = we assume 
that the qubit's state vector is |\l/(0))=a| ]) + b\ 1) with 
|a| 2 =|&| 2 =l/2, which in the Bloch vector language corre- 
sponds to the vector being in the xy plane. Experimen- 
tally this is achieved by initializing the qubit in one of 
the eigenstates of a z and applying a n/2 rotation about 
x or y axis at initial timeP 

In Free Induction Decay (FID) experiment we let the 
qubit evolve freely for time t, and then perform a mea- 
surement. Due to the noise in the Hamiltonian ([!]), the 
qubit state at the measurement time is 

|*(t)) = e -Mt/2 e -ifZfKt')dt' a \ ^ +e m/2 e i /*/J(f)«if 6 | |) 

so that the off-diagonal element of the qubit density ma- 
trix is (we use the units with K=l) 



(10) 



where p + -(0) — ab* . We quantify the qubit coherence 
using function W(t), defined in the following way: 



W(t) 



\(p+-(t))\ 

Kp+-(o))| 



(ii) 



We are now going to consider applying a certain num- 
ber n of ideal (<5-shaped) n pulses (about, for example, 



x axis) in the time interval if £ [0,i\. In the following, t 
will always denote the measurement time, and by W(t) 
we mean the coherence at time t with n pulses applied 
within this time. 

Since the n rotation about the x axis is given by 
exp(— z7T(T x /2) = — ia x , the qubit evolution operator with 
7r pulses applied at times ti,...,t n is 



U(t)=e-^'^ [n+p ^ ]dtf {-i& x )e 



^ [Q+/3(t')]dt' 

(12) 



The a x operators exchange the amplitudes of | f ) and 
I I) states of the qubit, and we arrive at the decoherence 
function under the action of the pulse sequence: 



W(t) 



cxpH / f3(t')f(t;t')dt') 
Jo 



(13) 



In this equation we have introduced the function f(t; t') 
which characterizes the pulse sequence: 



(14) 



fe=0 



where 0(t') is the Heaviside step function, to = and 
t n +x = t, the total evolution time. This function switches 
between 1 and —1 at the times at which the it pulses are 
applied and f(t; t') = for t' < and t' > t. In Fig. [T]d we 
show, as an example, a plot of f(t; if) for 2 pulse CPMG 
sequence. 

We denote the characteristic time of decay of W(t) as 
T2, defined by lnW / (T 2 ) = — 1. It depends on the pulse 
sequence applied during the qubit evolution, and this de- 
pendence on the number of pulses n and their spacing 
is the main subject of this paper. In most of the cases 
considered here the decoherence is not described by a 
simple exponential decay W(t)~ exp(— t/T 2 ). Such a de- 
cay law appears when the relevant dynamical time-scale 
of the environment (i.e. the noise autocorrelation time) 
is much smaller than T2 (the Markovian limit of qubit 
dynamics. This is not true for 1// noise which is corre- 
lated on a very long time-scale, and also for the RTN due 
to a slow fluctuator (with small 7). 



A. Pulse sequences aimed at suppressing the pure 
dephasing 

We denote the times at which the n pulses are applied 
by tk = 8kt with 0<i5fc<l, with k = l...n. The spacing 
of times tk in the sequences under consideration here is 
illustrated in Fig. [TJi. Application of a single ir pulse at 
t' = t/2 (n =1, S\ = 1/2) corresponds to the spin echo 
(SE) experiment. The tt pulse exchanges the amplitudes 
of the two states of the qubit, and the evolution during 
the remaining time period partially cancels the evolution 
before the pulse. More specifically, the echo sequence 
removes the influence of the noise frequencies lo smaller 
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FIG. 1: (a) The illustration of various pulse sequences with 
application times of n pulses marked. Spin echo (SE) is shown 
along with the PDD, CPMG, CDD and UDD sequences with 
n = 5 pulses (for CDD this corresponds to the 3 rd order of 
concatenation). 10 pulse UDD and CDD (4 th order of con- 
catenation) are also shown, (b) Function f(t;t') defined in 
Eq. (fl4j) for 2 pulse CPMG sequence. 



than 1/t (the quasi-static shifts of the qubit precession 
frequency). This, of course, gives the complete cancel- 
lation of the static randomness in qubit splittings f2 in 
the measurement of an ensemble of qubits (inhomoge- 
neous broadening) . In the considered here case of the SC 
qubits, one always deals with repeated measurements of 
a single qubit (a time ensemble) , and SE is a very efficient 
technique which suppresses the low-frequency part of the 
1/ / noise, leading to a su bstantia l increase in th e T% ti me 
in superconducting charge^^H anc i nux qubitsJ22l23 

The multiple-pulse extension of the echo is the Carr- 
Purcell-Meiboom-Gill (CPMG) sequence,^ defined by 
Sk=(k — l/2) jn. This sequence leads to periodic echo-like 
revivals of the coherence. While it was used for many 
years in NMR experiments performed on ensembles of 
spins, it was recently predicted that it should be highly 
effective at prolonging the coherence of a spin qubit in- 
teracting with the nuclear bathPSl Whereas CPMG is 
best understood as a refocusing sequence, in recent years 
a lot of attention has been devoted to the idea of dy- 
namic al decoup ling of the qubits from its environment by 
pulses LL i l 15 l 16 l 17 l l n this approach the pulses are used to 
average out the influence of the environment on the qubit, 
which can be fully achieved only in the limit of very fast 
repetition of pulses. Out of many types of dynamical 
decoupling (DD) protocols, we concentrate here on de- 
terministic periodic (PDD) sequence (see e.g. Ref. 33 for 



a comparison of more kinds of DD techniques applied to 
the spin bath problem). It is defined by 5k=k/(n + 1). 
Although it looks very similar to CPMG (the difference 
being only a small offset of the initial and final delay 
times), below we will show that CPMG visibly outper- 
forms PDD when considering the realistic small n. The 
key difference is that while in the limit of very fast ap- 
plication of pulses both PDD and CPMG decouple the 
qubit from the environment (during the whole time of 
the evolution), for realistic small n the CPMG sequence 
is much better at refocusing the coherence at the final 
time t' = t. 

Recently, a new family of DD protocols involving con- 
catenating (recursively embedding the sequences within 
themselves) has been proposedP^For the purpose of com- 
bating the pure dephasing, we will concentrate on con- 
catenations of the echo sequence, and for simplicity we 
will refer to it as simply CDD. The CDD sequence at Z-th 
order of concatenation and for total evolution time t is 
defined as CDD;(t). CDD (t) is free evolution for time 
t. The Z-th order of concatenation is then recursively 
defined by 

CDD,(t) = CDD,.! Qj -tt-CDD^! (0 -tt , (15) 

so that, CDDi(i) is the SE with a tt pulse at t/2 and 
CDD 2 (i) is the same as n — 2 CPMG sequence. For l>2 
the concatenations of the echo give us new sequences of 
non-trivially spaced tt pulses, the performance of which 
has been invest igated t heoretically in the case of the nu- 
clear spin bathpil221MI Note that the number of pulses 
for the Z-th order of concatenation is n ss 2 l . The CDD 
sequences were argued^ to be more tolerant to imple- 
mentation errors and more efficient (in terms of perfor- 
mance for the same number of pulses) than the PDD 
sequence. However, the theoretical comparisons^ with 
other DD protocols were done in the quantum mechani- 
cal setting, using the Magnus expansion of the evolution 
operators, or so-called "average Hamiltonian" theory.^ 
CDD is designed to cancel, with each order of concate- 
nation, successive orders of the qubit-bath interaction in 
the Magnus expansion. It also cancels successive orders 
of intra-bath interaction.^! However, it is not a priori 
clear whether the advantages of CDD are going to also 
hold for the case of dephasing due to classical noise. 

The most recent development in suppressing the pure 
depha sing was the introduction of a new sequence by 
Uhrig GsEsU which we term here UDD. This sequence was 
optimized for pure dephasing due to a bosonic environ- 
ment or classical Gaussian noise, but later its surprising 
universal character was discovered in a general quantum- 
mechanical setting.^ UDD is defined by 



S k = sin 2 [7r/c/(2n + 2)] 



(16) 



and in the next Section we will explain in what sense it is 
"optimal" for the case of the Gaussian noise. Originally 
the sequence was applied-^ to the case of the environment 
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characterized by an Ohmic spectral density of the noise 
having a sharp high-frequency cutoff, S(u>) ~ ljO(lj c — lj), 
with the ultra-violet cutoff uj c . Here we will analyze its 
performance for Gaussian noise with 1// spectral den- 
sity (with and without the ultra-violet cutoff in the spec- 
trum), and for classical non-Gaussian RTN. 



B. Realistic pulses 

We consider here the case of ideal, i.e. (5-shaped 7r 
pulses (so-called "bang-bang" or unbounded control). In 
reality, the pulses will have a finite duration t p and they 
might be imperfect, e.g. one can have pulse length or 
amplitude errors (leading to a wrong angle of rotation) 
or an off-resonance error, due to which the rotation oc- 
curs around a tilted axis. When these errors are system- 
atic, they can be suppressed by using composite pulses.^1 
Pulse shaping has also been used to counteract the ef- 
fect of the bath noise during the finite t p of the realistic 
pulse. Shapes of finite duration ir and tt/2 pulses were 
optimized to cancel the lowest order (in t v ) corrections 
due to interaction with arbitrary bath. 53 54 Optimization 
of control pulses was also considered for qubit coupled to 
a source of classical RTlNf^ or classical 1 // noisej^ and 
qubit interacting with a quantum two-level system.^ 

These works show that the realistic ir pulses can be 
made quite robust to both implementation errors and 
environmental noise, and treating them as 5-shaped is a 
good approximation as long as t p is larger than T m j n , the 
minimal interval between the pulses in a given sequence. 
For UDD this time scales with the number of pulses n as 
T min oc t/n 2 in contrast to t/n scaling for all the other se- 
quences under consideration, see Fig. [T^ for example with 
n = 10. In reality there is a lower limit on t p , related to 
the presence of higher energy levels in the full spectrum 
of the system. The inevitable higher-order pulse errors 
(unaccounted for by optimization) can also add up in a 
sequence with large n. Let us also mention that while 
CDD was showrP-3 to be robust against certain types of 
pulse errors by construction (the systematic errors being 
cancelled by successive concatenations), such an investi- 
gation has not been made in the case of UDD. ft is not 
known how sensitive is the performance of this sequence 
to, e.g., errors in timing of the pulses, which have to be 
spaced in a quite intricate fashion. These considerations 
lead us to concentrate on the case of rather small n. In- 
stead of looking at a "stroboscopic" limit of dynamical 
decoupling of the qubit by fast repetition of pulses, we 
start from the echo sequence and show how the decoher- 
ence changes as we increase n from one to ten. 



IV. GAUSSIAN NOISE. 



We write the decoherence function W(t) from Eq. ( 13 ) 



as 



defining the function x(U- I n the Gaussian approxi- 
mation, the average over noise can be performed using 
Eq. pj, and x(i) can be expressed through the spectral 
density of the noise S(lj) as 



x(t) 



du> 
2^ 



5(a,)|/(t;a,)| 2 = 



(18) 

where fit; lo) is the Fourier transform of f(t; t') with 
respect to f. The filter function F(wt) = ^\f{t;uj)\ 2 
encapsulate s the influence of the pulse sequence on 
decoherence! 24 * 58 ! In terms of times tk at which the pulses 
are applied (with to = and t n+ \ = t) we have 



k=Q 



(19) 



Analytical expressions for F(uit) for the sequences under 
consideration are given in Table |T] Let us note the exis- 
tence of the following sum rule for the filter functions: 



+oo 



du F(ut) 



f 2 (t;t')dt' = t 



(20) 



From this one can see that the pulse sequences cannot 
prolong the coherence time when the integral in Eq. ( 18 1 
is dominated by an initially flat S(lj) ~ 5(0) at low fre- 
quencies. We obtain then x(t) = 5(0)i/2 for all pulse se- 
quences at times t > with ujf being the frequency 
at which the noise spectrum starts to decay. 



The filter functions 



For free-induction decay (FID), as lj — > 0, F{u)t)/uu 2 — > 
t 2 /2. As a result, low frequency noise significantly con- 
tributes to x(u> an d for 1// noise with S(u) = Aq/u wc 
getP 



x(t) 



2ir 



In 



1 



(21) 



W(t) = e"*M 



(17) 



where u)\ t it the infra-red cutoff of the 1/f noise. In 
current experiments the microscopic cutoff is not reached, 
and u)- 1T is determined by the measurement procedure, 
i.e. LUij- ~ t m , where t m is the averaging time. Values of 
o;i r /27r~ 1 Hz have been reported for SC charge qubitsP^ 
This exposure to small-u> noise is already removed with 
the SE sequence, for which F(iot) ~ (wt) 4 for t <C 4/w. 
For 1/f noise this leads to a significant (by at least an 
order of magnitude) increase of th e observe d T 2 time in 
comparison to the FID experiment ! 5 I 11 I 2 Q 12U 

The PDD filter for lj < 2/t is F PDD (ujt) -cj 4 (lj 2 ) for 
odd (even) n (see Table [l} , so that only odd n sequence 
can suppress low-frequency noise. For larger frequencies, 
but smaller than 2n/t we have F PDD (ujt) - (ujt) 2 /(2n + 
2) 2 . On the other hand, the CPMG filter is propotional 
to lj 4 (lj g ) for odd (even) n, suppressing low frequency 



6 



Sequence 



FID 

SE 

PDD (odd n) 
CPMG (even n) 
CDD 

UDD 



m 

2 sin 2 f 
8 sin 4 § 

2tan2 2^+2 sifl2 f 

o ■ 4 z ■ 2 z / 2 z 

8 sin j^sin 3 /cos K 
2 2; + 1 sin 2 ^ r nLiSin 2 ^ 



1 V"^Tl 

2 | Z^fc = _. 



;fc + l 

lfc exp[f cos^] 



TABLE I: The expressions for filter functions for various pulse 
sequences. Here n is the number of pulses, and I is the order 
of concatenation for CDD (n~2 ). In the range of z < 2n, 
the filter function for UDD is very small, see Ref. 1181 In 
the formulas for even-n PDD and odd-n CPMG, sin 2 (z/2) is 
replaced with cos 2 (z/2). 



noise with S(u>) *~ l/w a with ev<2(4). Furthermore, for 
w < 2n/t we have F GPMG (ujt) ~ (ut) 4 / (2n) 4 . A small 
change of the initial and final interval between the pulses 
in comparison to PDD, leads to a more efficient high-pass 
filter of the noise. 

From the recursive definition of CDD we get in the Z-th 
order of concatenation 



/ (0 M = \ f {k ~ 1] 



(i 



iu>t/2 



), 



(22) 



with f(°}=i/u>(l — e lLJt ). From this the formula for the 
filter Fp BD (uit) given in Table [i] follows. For frequencies 
uj<4:/t we have 



F 



CDD 



21+2 



2('+!) 



(23) 



It is important to note that unlike in the case of the other 
pulse sequences, the frequency at which F(wt) becomes 
larger than 1 scales not as n/t, but as y/n/t (for large 
n). This is illustrated in Fig. [2j where the filters F(ujt) 
are shown for n= 10 for all the sequences under consid- 
eration. The CDD filter is the first to become large with 
increasing u>. The advantage of CDD over much simpler 
PDD and CPMG sequences is apparent only at very low 
frequencies - but it is the UDD filter which is clearly the 
best for w<c2n, see Fig. [2] 

The n pulse UDD sequence is derived by optimizing 
x(i), choosing 8k so that the first 2n + 1 terms in time 
expansion of x(t) around t=0 are zero (i.e. setting the 
first n terms in time expansi on of f(t; u) to zero). The 
physically meaningful solutiorPOH] to the resulting set of 
nonlinear equations is given by Eq. ( |16[ ). The resulting 
F(ujt) is 



F 



UDD 



(Wt) 



8(n+ l) 2 (iot\^+2 



112 



for u < 2 jt , (24) 



and F^ DD - [wt/{2n + 2)] 2 ™+ 2 < 1 for uj < (2n + 2) ft. 
Of all sequences considered here, the UDD gives the filter 
which most strongly suppresses the noise at low frequen- 
cies, as shown in Fig. [2] 




0.2 0.4 0.6 0.8 

co t/(2*n+2) 




0.4 0.6 0.8 

co t/(2*n+2) 



FIG. 2: (Color online) The filter functions F(ut) for SE and 
various pulse sequences with n — 10 pulses (for CDD it cor- 
responds to the 4th level of concatenation). The lower panel 
shows F(ut) in the logarithmic scale. 



Summarizing, the application of n pulses within time 
t effectively suppresses the noise power below frequency 
w n ~2n/t, with the UDD sequence being, by construction, 
the most efficient high-pass filter. Thus, at short time t 
or for large n only the high-frequency fluctuations with 
uj>2n/t contribute to x{t)- 



B. Gaussian 1// Q noise 

We concentrate now on the case of Gaussian 1// Q noise 
with spectral density S(w)=A] j +a /ui a , where 0.5 < a < 
1.5. The conditions under which the 1// noise origi- 
nating from multiple TLFs is Gaussian are discussed in 
Ref. [55] We first consider the case in which an ultra-violet 
cutoff uj c is present in the noise spectrum, as it was in- 
ferred for charge noise from experiments in Ref. |T1 . If we 
apply n pulses in time t such that 2n/t>uj c , all the noise 
is strongly suppressed, as shown see Fig. [3] Observation 
of an initially flat W(t) is a clear-cut signature of a finite 
cutoff. Therefore, pulse sequences can provide important 
insight into the noise spectrum. 

The decay of qubit coherence for various pulse se- 
quences is shown in Fig. [3] where we compare W(t) for 
various 5-pulse sequences with the echo (SE). The FID 
signal is not shown, since it depends on measurement- 
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FIG. 3: (Color online) The dependence of (a) W(t) and (b) 
x{t)——\nW(t) for SE and higher-order (n — 5) sequences for 
1/f noise with A /u c = 1. UDD and CPMG give W{t) « 1 
for £ < (2n + 2)/o> e - Echo signal (SE) (equivalent to PDD or 
CPMG with n—1) is also shown for comparison. 



specific infra-red cutoff lu- u -. However, in typical experi- 
mental situations M^ FID decays much faster than W SE . 
For a given n, PDD is clearly the least effective approach 
at all times. As expected, for short times f < 2n/w c , 
UDD is orders of magnitude better than the other pulse 
sequences, see Fig. [SJo. Thus, UDD is the ideal sequence 
for maintaining a low level of decoherence, i.e. high fi- 
delity, which is a necessary condition for quantum error 
correction. However, if the goal is simply to increase the 
characteristic decoherence time T 2 , defined by x(72)=l, 
then for a given n, the CPMG sequence is the best strat- 
egy- 
It is interesting to note that the CDD sequence does 
not offer strong advantages compared to CPMG and 
UDD. At short times it gives smaller x(t) compared 
to CPMG, but the difference is not as dramatic as in 
the case of UDD. At longer times, CPMG is better and 
gives larger T 2 . It seems that the benefits of using CDD 
expected in the regime of quantum bath dynamics are 
largely lost when dealing with classical noise. 

A different situation arises when the cutoff cannot be 
reached, i.e. we cannot apply pulses fast enough. Then, 



only the low-frequency part of Siui) can be suppressed, 
and the decay of W(t) is due to high-frequency tail of 
noise spectrum. For all the sequences under consider- 
ation it can be shown that to a very good approxima- 
tion (although slightly worse for PDD) the integration of 
Eq. (fT8|) gives 



x(t) 



C a (Apt) 1 
2ir n a 



(25) 



Here C a is a sequence-specific constant of the order of 
one. This result can be quickly established for CPMG in 
the following way. Its filter function F(z = ivt) consist of 
a periodic train of narrow peaks, with periods z p = 2im. 
We approximate these peaks by square steps of width 
Az, which can be derived from the sum rule for F(z) 
(Eq. (20 1). Then, to the lowest order in small quantity 
Az/ z p = l/n we get for that for CPMG we have 



7T 2 (27T 



1 



k=l 



(k — 1/2) c 



(26) 



which for a = 1 gives Ci = 7C(3)/tt 2 ~ 0.85 (where ((z) 
is the Ricmann zeta function). On the other hand, for 
UDD one can use a seemingly crude approximation of 
F,y DD (z) = a n 6(z-(2n + 2)), with a constant a n . This is 
motivated by the fact that the UDD filter for z > 2n + 2 
looks similar to a random signal, since it is given by a 
sum of periodic functions with non-commensurate pe- 



riods (see Eqs. (16) and (19l). From Eq. (201 we get 
o-n = 7t(n + 1), and this value gives a good agreement 
with results of numerical integration using the exact form 
of the filter. The analytical approximation gives then 
C a = Ti 2 /[(«+ 1)2 Q+1 ] and n should be replaced by n+ 1 
in the denominator in Eq. (25 1. From these formulas, 



and from numerical calculations confirming their accu- 
racy, we find that CPMG marginally outperforms the 
other sequences (more visibly for larger values of a), so 
it is enough to implement this simple sequence to prolong 
qubit coherence in this regime. 



C. Noise spectroscopy using the pulse sequences 

The time dependence of qubit coherence under exter- 
nal pulse sequences can be used as a spectroscopic tool 
for extracting the noise spectrum contributing to dephas- 
ing. The idea of using qubit energy relaxation for noise 
spectroscopy was introduced in Ref. 1401 and it has since 
been realized experimentally. 41 Here we propose a quan- 
titative method for extracting the moments of the noise 
contributing to pure dephasing, which can be different 
than the noise leading to the energy relaxation. In par- 

2n+2 



ticular, for t<l/w c we have x (t) ~ t 2n+2 M 2n (see 
Eq. ( 24 1 ) , where Mj~ = J oj k S(u>)du> is the k-th moment 



of the spectral density. From the moments Mj~ one can, 
in principle, reconstruct the noise spectrum. For SE and 
two-pulse UDD (equivalent to two-pulse CPMG) we get, 
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respectively, X ~ M 2 t A /32n and X « M 4 t 6 /1024tt. The 
observation of exp(— i 4 ) and exp(— t 6 ) decays of W(t) for 
these one and two pulse sequences will be a signature of 
the presence of finite uj c . 

Fulfilling the condition 2n/t > uj c might however be 
experimentally challenging. As we discussed in Section 



III B in reality the pulse time t p has a lower bound, and 
n can also be limited by accumulation of errors in a long 
and complicated sequence. Assuming r p is the limiting 
factor, if 2/t p ^>u} c , it is possible for the filter function to 
"reach the cutoff" , and then the previous considerations 
hold. 



V. NON-GAUSSIAN RANDOM TELEGRAPH 
NOISE. 

The comparison between the experiment^ and theory^ 
clearly shows that in charge qubits the decoherence can 
be dominated by coupling to a single classical fluctua- 
tor, which is a source of the RTN ((3(t) = v£(t) with 
£(t) switching between ±1/2 wi th rat e 7). Two regimes 
of decoherence can be identified ] 6 ! 9 ! 29 ! the strong (weak) 
coupling regime in which g=v/j^>l (g<^l). For g <C 1 
we are in the "motional narrowing" regime: the fluctu- 
ator is switching so fast that its influence on the qubit 
is averaging itself out, leading to large T 2 . Furthermore, 
since on the relevant time-scale the qubit receives a large 
number of "phase kicks" from the fluctuator (with typ- 
ical size of v/j), the effective noise affecting the pure 
dephasing dynamics is approximately Gaussian. On the 
other hand, for g 3> 1 one expects short decoherence time 
with str ongly non-Gaussian features in time-dependence 

of wit)0xto 

We have studied the effect of the pulse sequences on 
qubit decoherence using both numerical simulations of 
the RTN and the Gaussian approximation, in which we 
plug the Lorentizan first spectral density of the RTN, 
Eq. ([8), into Eq. Q. The results for W(t) in both 
coupling regimes are shown in Fig. [4} For <?<Cl, the effect 
of pulses is marginal, i.e. one has to apply a large number 
of pulses to obtain a visible effect. On the other hand, 
in the strong coupling regime, application of even a few 
pulses substantially increases the coherence time. Similar 
to the case of the Gaussian 1 // Q noise without a cutoff, 
the CPMG sequence is the better practical approach. 

For g^>l there are strong deviations from Gaussian 
behavior in the SE signal (see Fig.|4^), and the shape of 
W(t) containing the characteristic plat eaus h as been de- 
rived using various analytical methods! 7 ! 9 ! 61 ! The values 
of v and 7 can be inferred from the position and height 
of the first plateauP However, one can see in Fig. [4^ that 
as we apply more pulses the deviation between the sim- 
ulation of the exact RTN and the Gaussian approxima- 
tion decreases. Therefore, with increasing n the simple 
analytical results following from Eq. ( 18 1 become more 
accurate, i.e. the non-Gaussian effects are suppressed by 
pulses. 
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FIG. 4: (Color online) The decay of qubit coherence W(t) for 
a single TLF coupled to the qubit for (a) strong coupling g = 
10 and (b) weak coupling g — 0.2. The results of the simulation 
of the RTN are shown as symbols, and the calculations in 
Gaussian (G) approximation are shown as lines. CDD4 (with 
n = 10) gives practically the same result as UDD, and thus 
it is not shown. For g = 0.2 the Gaussian approximation 
agrees very well with the exact results. For clarity, only SE 
and CPMG with n=10 are shown in (b). The FID signal 
for strong coupling (not shown) is an oscillating function for 
which the SE signal is an envelope. 



The explanation of the improvement of Gaussian ap- 
proximation with increasing n in the strong coupling 
regime is the following: the deviation between the exact 
result for RTN and Gaussian approximation arises from 
higher order noise correlators in the cumulant expansion 
of W(t): 



\nW{t) = -x [2) {t)- X ^{t) + 



(4)/ 



(27) 



with ~<?"i an d X^ given by the expression following 
from the Gaussian approximation, Eq. (18). The ratio 
R(t)= X^\t)/x (t) can be used as a measure of the 
importance of non-Gaussian effects. We have calculated 
it for various pulse sequences, finding that x' 4 ^ is more 
strongly suppressed by pulses than x^ 2 \ so that while 
the coherence time T 2 is extended with increasing n, the 
time-scale on which the non-Gaussian effects are negli- 
gible grows even faster. The details of the calculations 
are given in Appendix [A] In Fig. [5 we show that R(t) 
remains small for a longer time with the application of 
more pulses. The CPMG sequence is better than UDD 
at suppressing R(t), which should not be surprising in 
the light of the fact that UDD is optimized to make only 



X 



(2) 



as small as possible. Evidently it suppresses x^ less 
efficiently than the CPMG sequence. 

In the Gaussian approach, using the analytical approx- 
imations outlined in Section llVBI we find that for CPMG 
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FIG. 5: (Color online) The ratio R(t) = X (4) /x (2) plotted for 
sequences having n — 0...5 pulses, with RTN characterized by 
7=1 and energy of interaction with the qubit v = 10 (so that 
g = 10). The solid lines correspond (from left to right) to 
FID, SE, and CPMG with n = 3,4, 5. The dashes lines corre- 
spond to the UDD sequence with n = 3, 4, 5. With increasing 
number of pulses, the time at which the non-Gaussian effects 
start to be quantitatively important (when R(t) > 1) becomes 
larger, eventually surpassing the T2 time of the qubit. 



VI. CONCLUSIONS 

We have analyzed the influence on various pulse se- 
quences on pure dephasing of a qubit affected by clas- 
sical noise, with emphasis on types of noise relevant for 
superconducting qubits. We have shown that successive 
higher order pulse sequences lead to an improvement of 
coherence time for both Gaussian l/f a noise and RTN. 
We have found that in the presence of a hard ultra-violet 
cutoff in the Gaussian noise spectrum, the UDD sequence 
is optimal for suppressing initial decoherence. However, 
if one can not "reach the cutoff" , the CPMG sequence 
is the best practical approach. This is also true in the 
case of a single fluctuator coupled to the qubit. There, 
the application of large n-sequences decreases the devia- 
tion between exact (non-Gaussian) theory and Gaussian 
approximation. For both 1/f noise and RTN we pre- 
dict substantial practical enhancement in SC qubit co- 
herence under the CPMG pulse sequence. Furthermore, 
a detailed experimental investigation of the noise mech- 
anisms operational in different SC circuits and samples 
becomes possible using the UDD pulse sequence, which 
allows one to gather quantitative information about low 
frequency noise contributing to dephasing. 



in the large n regime we have 



x(t)< 



24 n 2 ' 



(28) 



-{"ft — n), jt^>n, 



and for UDD the numerical coefficient in the first equa- 
tion is larger by a factor of about 1.5 and n should be 
replaced by n + 1 in the denominator. The first formula 
holds when the filter function F(u>t) suppresses the low- 
frequency (u < 27) flat part of the Lorentzian spectral 
density, and only the jv 2 /lu 2 tail contributes to If 
g^> 1, most of the decoherence occurs for x(t) ^t 3 , an d 
T2 falls within this regime. Then the application of more 
pulses is effective as it decreases the coefficient of t 3 , re- 
sulting in 



T 2 



2ti 2 / 3 



19 



2/3 



ff » 1 



(29) 



On the other hand, for g <C 1, most of decoherence occurs 
in the long time (-ft 3> n), x(i)~i regime, in which the 
largest contribution to \ comes from the flat part of the 
spectral density. Then, adding a few pulses only prolongs 
the initial short-time x{t) ~ t 3 behavior, with marginal 
effect on the decoherence time (T2 « 8/-fg 2 ). However, 
the initial decoherence (-ft <C n) is suppressed as before, 
resulting in improvement of fidelity at short times. The 
pulses affect the T 2 time only when we apply n > 8/g 2 
pulses, extending the t 3 regime so that T2 falls within it. 
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APPENDIX A: CALCULATION OF THE 
FOURTH CUMULANT OF THE RANDOM 
TELEGRAPH NOISE 



We define the phase $(t): 



*(*)=« / i{t')f{t; t')dt' , 
Jo 

so that the decoherence funtion W(t) is given by 



W(t) = (e" l * (t) ) = exp 



(E 



H)" 



k=l 



(Al) 



(A2) 



where we have written it using the cumulant expansion. 62 
The cumulants Ck vanish for k > 2 if the statistics of Q(t) 
is Gaussian. They can be written in terms of moments 
Mj~(t) — (Q k (t)), with the first two non- vanishing ones 
(we assume (£) = 0, so that the odd-fc moments and cu- 
mulants vanish) given by 



C 2 (t) = M*(t) , 

Ci(t) = Mf(t)-3[Mi{t)} 2 



(A3) 
(A4) 
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We also define 



lnW(t) = x(t)=X> (2fc) W 



(A5) 



fc=i 



In the Gaussian approximation the only non-zero term 
in the above expansion is x^ — For the RTN 

the higher-order terms do not vanish, e.g. we have 

In order to calculate x (t) we need to understand the 
structure of the higher-order correlation functions of 
Following Ref. [7] we write (3{t) as 



m = < (-i) n(o ^ , 



(A6) 



with the result being the same for h > h (since n/h, h) — 
n(ti, h)), so that we recover Eq. Q. In an analogous way 
we can calculate the four-point correlation function, but 
now the ordering of time arguments will matter. Assum- 
ing h > t% > h > £4 we get 



16 



. e -27(*i-*a)g-27(*3-f4) 



(A9) 

Time ordering is crucial here. For any other ordering 
we have to permute the times on the right-hand side. 
However, we deal here with multiple integrals of the form 



where £o=£(0)= ±1/2 is the initial condition, and 
n(ii,i2) is the random variable giving us the number of 
flips between times t\ and h- From definition of the RTN 
process we have 



((— l)"(*!' t2 )) = g — 2 Tl*l— *2| 



(A7) 



The two-point correlation function of the noise can be 
written for t\ > h as 

(P(h)(3(h)) = «^a (( _ 1) n(o,* 1 )(_i)»(o,fc) ) 
..2 



V _/(_- i \n(t 2 ,t 1 )\ _ ^ e -27|t!-t 2 | 

4 4 



(A8) 



( / drar)f(T)Y 
Jo 

dti... / dt k {Z(t 1 )..4(t k ))f(t 1 )...f(t k ) ,(A10) 
Jo 



where we have used the simplified notation f{h) = 
f(t;ti). The integration region (the fc-cube) can be di- 
vided into k\ simplexes, each with a definite ordering rela- 
tion between all the times. The integration variables can 
be relabeled in each integration region, and we obtain 



( / drC(r)/(r))' 



k\ / dh 



dto 



tk-1 



dt k (ah)..4(tk)}f(h)...f(t n ) 



With this formula we get for the moments: 



Mi = ~ 



dh / dt 2 e- 27( * 1 - t2) /(£i)/fe) 



dh / dt 2 / dh / dhe- 2 ^- u - +t ^f(h)f(h)f(h)f(h) 



(All) 

(A12) 
(A13) 



where the formula for Mf — 2x^ is simply a different way of obtaining the Gaussian result from Eq. (18 1 with 
Lorentzian spectral density. Using Eq. (A4| we obtain the fourth cumulant C^(t). For small number of pulses n, e.g. 
for FID, SE, and n = 2 CPMG/UDD (labeled hereafter as CP2) we get 

(A14) 

-T t (l + 2 7 i)-13 X ) , (A15) 
. . 24e- 57t/2 - 12e- 27t + 24 e - 37 ' /2 + 8e" 7t + 8e" 7 * /2 — 21^ . (A16) 



C 4 (t) FID = 


3<? 4 /„ 


e" 47 * 


Cf(t) = 




e" 47 * 


C 4 CP2 (i) = 




+ 2e" 




+e -47t _ 8e - 


-7 7 t/2 



24e 



-37* 



-37t/2 



- 7 t/2 



The analytical expressions for larger n become cumbe- 



some, and we resort to numerical evaluation of C4. The 
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results for the ratio of the cumulants 



R(t) 



x (4) W 



1 C A {t) 



x W(t) l2C a (*) 



(A17) 



up to n = 5 are presented in Fig. [5] 

The fact that with increasing n the higher order cu- 
mulants are suppressed more strongly than the Gaussian 
C2(t) can be understood in the following way. ~x 2k is 



proportional to 2k— fold time integral of a noise corre- 
lation function multiplied by 2k functions f{U), each of 
them alternating between ±1. Under the multiple in- 
tegral, and for large n, the sign of the product of /(t») 
switches multiple times, and with increasing order 2k the 
whole expression is effectively averaged out by the filter 
functions. 
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